Field and current distributions and ac losses in a bifilar stack of superconducting strips 
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In this paper I first analytically calculate the magnetic-held and sheet-current distributions gen- 
erated in an infinite stack of thin superconducting strips of thickness d, width 2a 3> d, and arbitrary 
separation D when adjacent strips carry net current of magnitude I in opposite directions. Each 
strip is assumed to have uniform critical current density J c , critical sheet-current density K c — J c d, 
and critical current I c = 2aK c , and the distribution of the current density within each strip is as- 
sumed to obey critical-state theory. I then derive expressions for the ac losses due to magnetic-flux 
penetration both from the strip edges and from the top and bottom of each strip, and I express the 
results in terms of integrals involving the perpendicular and parallel components of the magnetic 
field. After numerically evaluating the ac losses for typical dimensions, I present analytic expressions 
from which the losses can be estimated. 

PACS numbers: 74.25.Sv,74.78.Bz, 74.25. Op,74.25.Nf 



o 
o 

u 

C3 



I 

C 

o 
o 



> 

00 
00 

o 
o 

oo 
o 



X 



I. INTRODUCTION 

The magnetic-field and sheet-current distributions 
generated in an infinite stack of superconducting strips, 
all carrying current in the same direction, were calcu- 
lated analytically in Refs. 0] and using an extension 
of a method first introduced by Mawatari^ The re- 
sults were then put to use to calculate the hysteretic ac 
losses. Such calculations can be applied to estimate the 
ac losses in pancake coils wound from long lengths of 
second-generation (2G) high-temperature superconduct- 
ing tapes4£ 

Recently noninductive coils with bifilar windings (in 
which adjacent tapes carry current in opposite direc- 
tions) have been fabricated using 2G superconducting 
tapes for use in superconduct ing fault-current limiters 
(SCFCLs)&I£&i° References [llf and 01 give excellent 
reviews of this topic. Each tape consists of a supercon- 
ducting film of the order of 1 /im in thickness and 1 cm 
in width. A thin insulating buffer layer separates the 
film from the underlying metallic base, which is typically 
tens of /im in thickness. Usually surrounding this struc- 
ture is a normal- metal (e.g., copper) stabilizer, such that 
the total tape thickness is a fraction of 1 mm. Account- 
ing for the thickness of the insulation between tapes, the 
spacing D between the superconducting films in adjacent 
tapes is of the order of a few mm. Since the interleaving 
tapes in such coils carry current in opposite directions, 
the current-generated magnetic fields arc localized within 
the windings and decay very rapidly outside the coil. 

To determine the current distribution within a bifilar 
stack of superconducting strips is not trivial. Roughly 
speaking, when a strip carries an increasing current, the 
current density just slightly exceeds the the critical cur- 
rent density J c only at the strip edges, where vortices 
penetrate and carry perpendicular magnetic flux through 
the strip. The middle portions carry a current density 
less than J c , and therefore no perpendicular magnetic 
flux penetrates the strip there. However, there is a mag- 



netic field parallel to the strip, and vortices can enter the 
top and bottom of the strip, remaining nearly parallel to 
the surfaces. 

In Sec. II, I present solutions for the current-density 
and magnetic-field distributions in an infinite stack of 
superconducting strips, each strip carrying a current of 
magnitude I but adjacent strips carrying currents in op- 
posite directions. In Sec. Ill, I show how to calculate 
the hysteretic ac losses generated by magnetic flux pen- 
etrating both from the strip edges and from the top and 
bottom surfaces. I give a brief summary and discuss the 
results in Sec. IV. 



II. INFINITE BIFILAR STACK 

As in many earlier calculations of the properties of 
thin-film superconductors, I consider only high-K type- 
II superconductors and assume for simplicity that the 
magnitude of the self-field H at the film edges or top 
and bottom surfaces is typically much larger than H c \, 
such that the magnetic induction in the superconducting 
film is given to good approximation by B = [1qH . For 
small currents this assumption leads to an overestimate 
of both the degree of magnetic flux penetration and the 
corresponding ac losses. I also will treat the quasi-static 
penetration of vortices into the film using critical-state 
theory^ parameterized by a critical current density J c 
that is independent of the local magnetic induction. 

Figure 1 shows the film geometry under consideration, 
an infinite stack of superconducting strips of width 2a, 
thickness d -c 2a, and infinite length parallel to the z 
axis, equally spaced along the y axis with separation D. 
The strips are assumed to be identical, characterized by a 
uniform critical current density J c , critical sheet-current 
density K c — dJ c , and critical current I c — 2adJ c — 
2aK c . 

Since the magnetic-field and current-density distribu- 
tions depend only upon the coordinates x and y it is con- 
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FIG. 1: Infinite bifilar stack: a stack of thin (d -C 2a) su- 
perconducting strips of infinite length in the z direction, with 
those at y = 0, ±2D, ±4D, ... carrying current / in the +z 
direction and those at y = ±D,±3D, ... carrying the same 
current in the —z direction. 



vcnicnt to describe the magnetic field outside the strips 
as H(C) = H y {x, y)+iH x {x 1 y), which is an analytic func- 
tion of the complex variable Q = x + iy. The Biot-Savart 
law for the complex fiel d 14 ' 15 i 16 can be expressed as 



1 f +a 

H{0 = rJ_„ 



du 



K zn {u) 



c 



inD 



(1) 



where K zn {x) is the sheet-current density in layer n cen- 
tered at (x,y) = (0,nD). The currents in a bifilar stack 
are distributed such that for n = 0, ±2, ±4, ±6, the 
sheet-current density is K zn (u) = K z {u), the same as in 
the layer n = 0, and for n = ±1, ±3, ±5, the sheet- 
current density is K zn {u) = —K z {u). The resulting sum 
in Eq. ([1]) can be evaluated, yielding 



1 f +a 

no = w j_ a d u 



K z {u) 



2D ' ''" 



sinh[7r(C - u)/D] 
K z {u) sinh(7rC/£>) cosh(7ra/£>) 
sinh 2 (7rC/D) -smh 2 (nu/D) 



(2) 
• (3) 



where the second expression follows from the symmetry 
that K z {— u) = K z (u). Note that the complex field has 
the desired properties that H(( + iD) — —H(() [i.e., 
H x (x,y+D) = -H x (x 7 y) and H v (x, y+D) = -H y (x,y)] 
and U[H{x + iD/2) = H v (x, D/2) = 0. 

We seek the solution for which, when a current / is 
first applied to the strips in the stack, the sheet-current 
density K z is equal to K c within bands of width (a — c) 
at the edges but obeys K z < K c in the middle region, 
|x| < c. The simplest way to obtain this solution is 
to change variables as follows: ( — (D/ir) sinh(7rC/D), 



u = (D/n) sinh(7ru/Z?), a — (D/n) sinh(7ra/D), and 
c = {D/ir) sinh(7rc/D), similar to the procedure used by 
Mawatari in Refs. H [lj andlH. Using K z (u) = K z (u) 
and H(C) = W(C), this yields 



27T La C ~ U 2 



(4) 



This is the Biot-Savart law for an isolated strip carrying a 
current density K z {u) = K z (—u), for which the solution 
is known to b o 19 i 20 ' 21 



Kr 



■ tanh 



~2 ' 



(5) 



where here and in later similar expressions \J C 2 — c 2 is 

shorthand for (£ — c) 1 / 2 (C + c) 1 / 2 . Thus the desired com- 
plex field is 



7T V sinh 2 (7rC/£>) -sinh 2 (7rc/D) V 



The corresponding complex potential, 

liD/2 



g(o = [ c neve, 

JiD/2 



(7) 



can be evaluated numerically. Contours of constant 
$tG(x + iy) (see Fig. 2) correspond to magnetic field lines. 

Taking the real and imaginary parts of Tl(x ± ie) in 
Eq. © and using K z {x) — H x {x — ie) — H x (x + ie) = 
2H X (x — ie) yields 



H v {x,0) 




(8) 



sinh 2 (7ra;/-D) — sinh 2 (7rc/D) 



m 1/ sinh 2 (7ra/D) - sinh 2 (7rc/D) ' 

c < \x\ < a, (9) 



— - tanh 1 



/sinh 2 (7ra/D) - sinh 2 (vrc/D) 
sinh (ttx/D) — sinh (ttc/D) 



\x\ > a, 



(10) 



K z {x) 



2K r 



tan 



/sinh 2 (7ra/L>) - sinh 2 (ttc/D) 
sinh {ttc/D) — sinh {ttx/D) 



\x\ < c, 
K Cl c < \x\ < a. 



(11) 
(12) 



The relationship between c/a and the current / carried 
by one of the strips is determined by the integral / = 
K z {x)dx, which can be expressed as 



I 2 

I r Tra 



, smh z {TTc/D) -smh. z {ttx/D) , 
tan" 1 \ I — -V; — -gi — '—^dx, 



sinh {na/D) — sinh {ttc/D) 



(13) 
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FIG. 2: Contours of constant $lQ(x + iy) calculated from Eqs. 
([6]) and (0 for a stack of thin superconducting strips of width 
2a (thick lines) and spacing D. The contours correspond to 
magnetic field lines, and in this figure when the current in 
is the +z direction, the field lines circulate in a counterclock- 
wise direction around the strips centered at (x, y) = (0, —2D), 
(0,0), and (0,2D), and in a clockwise direction around the 
strips centered at (0, —D) and (0, D). Here, D = a = 2c. 
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FIG. 3: Real (solid) and imaginary (dashed) parts of $tH(x — 
ie) — H y (x,0) + iH x (x, — e) calculated from Eq. © or Eqs. 
(l8"))- (|12p for a stack of thin superconducting strips of width 2a 
and spacing D. Here, o = 1, c = 0.5, and D = 1. 



where I c — 2aK c . Shown in Fig. 4 are plots of c/a vs 
{I/Icf for D/a = 0.1, 0.3, 1, 3, and 10. For D/a = 10, 
the plot of c/a vs (I/I c ) 2 is nearly indistinguishable from 
c/a = \J\ — (I/I c ) 2 , the result known for an isolated 
strip i 19 ' 20 ' 21 On the other hand, for small values of D/a, 
the plot of c/a vs (I/I c ) 2 can be calculated to good ap- 
proximation by 

sinh(7rc/L>) = sinh(7ra/L>) cos[(7r/2)(/// c )], (14) 




FIG. 4: Plots of c/a vs {I/I c ) 2 (solid) determined from Eq. 
p3|) for D/a = 0.1, 0.3, 1, 3, and 10 (top to bottom). 
The dashed lines show the corresponding linear slopes for 
(I/Q 2 < 1 [Eq. (IS)]. 
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FIG. 5: Plots of K z (x)/K c vs x/a for I/I c = 0.5 determined 
from Eqs. CLE]), dHJ), and ([T3j) for D/a = 0.1 (c/a = 0.989), 1 
(c/a = 0.926), and 10 (c/a = 0.868). 



such that the value of c/a is very close to 1 except when 
I/I c is very close to 1. This behavior occurs because the 
current density K z (x) in the region |x| < c is practically 
constant and nearly equal to I /2a for a wide range of 
subcritical values of /, as shown by the plot of K z (x) vs 
x for D/a = 0.1 in Fig. 5. The curve of c/a vs (I / I c ) 2 for 
D/a = 0.1 is indistinguishable from that obtained from 
Eq. fli]). For all values of D/a, S = (a - c) = aa(I/I c ) 2 
for small values of (J/J c ), where Eq. (fTS"]) yields 



i/D 



2tanh(7ra/L>) lK(k) 



tt/2 



(15) 



and K(k) is the complete elliptic integral of the first kind 
of modulus k — tanh(7ra/D). This behavior is illustrated 
by the dashed lines in Fig. 4. In the limit as D/a — > oo, 
Eq. (IS]) yields a = 1/2, and for D/a « 1, a = TrD/8a. 
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III. AC LOSSES 

A secondary goal in this paper is to calculate the 
hysteretic ac losses in a bifilar stack of superconduct- 
ing strips at a frequency / = 1/T that is sufficiently 
low that eddy-current losses are negligible and the losses 
can be calculated using a quasistatic approach^ The 
solutions for H(x,y) derived in Sec. II can be used to 
calculate Q', the energy dissipated per cycle per unit 
length in each strip. Consider time t = 0, when the cur- 
rent / has its maximum amplitude in the z direction, 
the magnetic- field distribution is given by H(x, y) — 
xH x (x,y) + yHy(x,y), and the magnetic induction is 
B(x,y,0) = fioH(x,y). Half a cycle earlier, at time 
t = — T/2, when the current is in the opposite direction, 
B(x, y, —T/2) = —)1qH{x, y). The loss per cycle per unit 
length Q 1 is twice the loss in the half cycle —T/2 < t < 0. 
Thus 

/0 pa pd/2 

dt dx dyJ z {x,y,t)E z {x,y,t). (16) 
-T/2 J -a J-d/2 

According to critical-state theory^ during this time in- 
terval, E z is nonzero only where J z is just above J c , such 
that J z can be replaced by J c in Eq. ([To]) , but the integral 
is to be carried out only over those portions of the cross 
section where E z (x,y,t) > 0. Note that E z (0,0,t) = 
when the current amplitude / is less than I c . Next, let us 
use Faraday's law in the form <f dl ■ E — — J dS ■ dB/dt, 
where the surface S consists of two rectangular parts of 
length L z , one with width x extending from the origin to 
(x,0) and the other of width y extending from (x,0) to 
(x, y). Integration of Faraday's law thus yields 

(17) 

Substituting this expression into Eq. (flB|) . integrating 
over time, noting that B(x,y,0) — B(x, y, — T/2) — 
2/ioH(x 1 y), and making use of the symmetry that the 
losses in the left and right halves of the strip are the 
same, we obtain 

Q' = 8n J c dx dy dx'H y (x',0) 

JO J-d/2 JO 

pa pd/2 py 

-8/x J c l dx I dy I dy'H x (x,y') 

Jo J-d/2 Jo 

= Q'e+Q'tb- (18) 

From this expression we see that there are in general 
two important contributions to Q' in superconducting 
strips. The first of these is Q' e , the dissipation due to 
magnetic flux in the form of vortex or antivortex seg- 
ments transporting flux density noH y in from the edges 
of the strip. The second contribution is Q' tb , the dis- 
sipation due to magnetic flux in the form of vortex or 
antivortex segments transporting flux density /ioH x in 
from the top and bottom surfaces of the strip. 




I/Ic 

FIG. 6: Plots of the dimensionless function describing edge 
losses L e vs F = I/I c for D/a — 0.1, 0.3, 1, and 10 (bottom 
to top). The solid curves were obtained from Eq. (|2ip and 
the dashed lines from Eqs. (|15[1 and (|26[1 . On this scale the 
solid curve for D/a — 10 is indistinguishable from the Norris 
result 1/2 [Eq. ([22]) ] for an isolated strip. 19 



A. Edge losses 

The edge losses Q' e are most conveniently calculated 
from the fundamental equatioi>^£ i 18 ' 19 i 22 i 23 i 24 i 25 

Q' e = 8fM K c [ (a-x)H y (x,0), (19) 

J c 

which is obtained from the first term in Eq. (|18j) by par- 
tial integration. Using the expression for H y (x,0) from 
Eq. ([9]), we obtain 

Q' e = ^I 2 c L e {F), (20) 
where F — I / I c and the dimensionless function L e is 

2 f a , , / sin h 2 2l3L_sinh 2 ^ 

L e = — - / (a — xjtanh \ Tr-^-dx, 

na 2 J c y sinh 2 £f - sinh 2 ^ 

(21) 

which is plotted as the solid curves in Fig. 6 for D/a = 

0. 1, 0.3, 1, and 10. 

In the limit as D/a — > oo, the results correspond to 
those for an isolated strip, and L e (F) reduces to the func- 
tion L 2 (F) derived by Norris^ 

L ^ F)= (1-F)\n(l-F) + (1 + F)ln(l+F)-F 2 

IT 

The maximum edge losses occur when F = 1 or I = I c . 
These can be calculated by setting c = in Eq. (|2"Tj) . In 
the limit as D/a — > oo, this equation yields 

L e (l) = (21n2-l)/7r = 0.123, (23) 

as obtained by Norris for an isolated stripJ^ For D/a <c 

1, on the other hand, Eq. (|2"Tj) yields 
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The smallness of this result arises because H y (x,0) f» 
(Kefir) tanh - {exp[— ir(a — x)/D]} when ira/D 3> 1, 
such that, except for the logarithmic divergence at 
x = a, H y (x,0) is exponentially small, i.e., H y (x,0) « 
(K c /ir) exp[— ir(a — x)/D], over most of the range of in- 
tegration in Eq. (12"T]) . The following interpolation func- 
tion between large and small values of D /a approximates 
L e (l) with a maximum error of 5%: 



Mi) 



0.123 



[l + 2.1(a/L>) 5 / 2 ] 4 / 5 ' 



(25) 



Expanding the right-hand side of Eq. (|2"Tj) to lowest 
order in powers of F = I / I c yields the approximation 



2 r 4 



L e (F) « -a'F 



(26) 



where a is given by Eq. (|15p . For large values of -D/a, 
this equation yields 



L e (F) » F 4 /6^, 



(27) 



as found by Norris,— and for small values of D/a, Eq. 
(H3) yields 



v ; 96 V a/ 



(28) 



Note from the solid curves and dashed lines in Fig. 6 that 
the approximation given in Eqs. (j!5[) and (|26[) provides a 
good estimate of the edge losses over a remarkably large 
range of values of F = I/I c - 



B. Top-and-bottom losses 

In the theoretical analysis of the ac losses in thin 
films ) 19 ' 20 usually only the losses due to vortex and an- 
tivortex motion in from the edges are taken into ac- 
count. These losses, represented by the term Q' e , dom- 
inate in isolated films when the ac current amplitude / 
approaches I c , for then the entering vortices travel an ap- 
preciable fraction of the strip width 2a during each cycle. 
See Eq. (fl9|) . On the other hand, vortices and antivor- 
tices entering from the top and bottom of the strips can 
travel at most a distance d/2, and when d <C a, it makes 
sense to ignore the top-and-bottom losses, represented 
by the term Q' tb . However, as seen in the above section, 
the edge losses are proportional to F A and also are much 
reduced for small values of D/a. Since for small F the 
top-and-bottom losses vary as a lower power of F, it is 
important to determine the conditions under which these 
losses exceed the edge losses. 

To evaluate Q' tb , we make use of the fact that, although 
in general J z (x,y) = dH y (x,y)/dx - dH x (x,y)/dy, in 
thin films the second term is far larger in magnitude. 
Therefore, to excellent approximation, 

H x {x,y) = - J c (y - y p ), y p < y < d/2, (29) 
= 0, -yp < y < y p , (30) 
= -J c (y + y P ),-d/2 < y < -y p , (31) 



where H x (x,-y) = -H x (x,y). Here y p (x) = d/2 - 
H x (x, -d/2)/J c for |x| < c, where H x (x, -d/2) < J c d/2, 
or y p (x) = for c < |x| < a, where H x (x, —d/2) = J c d/2. 
Carrying out the second integral in Eq. (|18[) . we obtain 



Qtb ~ sT c 



\H x (x,±d/2)\ 3 dx. 



(32) 



This result is expected, because when a type-II supercon- 
ductor is subjected to a parallel ac field of amplitude Hq, 
the hysteretic ac loss per unit area per cycle is known to 
be±2 



Qa = 



2/ioff 3 
3J C 



(33) 



The result for Q' tb in Eq. ([32]) corresponds to replacing Hq 
by \H x (x,±d/2)\ and integrating Eq. (f3"3"|) over the top 
and bottom of the film. Since we are here considering the 
case that d <C a, we can evaluate Eq. (|3"2"]) by replacing 
\H x (x,±d/2)\ by H x (x,-e) = K z {x)/2, where K z {x) is 
given by Eqs. (jTTJ) and (TT2"|) . The top and bottom losses 
therefore can be expressed as 



Q'tb = IMiI c Ltb{F), 



(34) 



where F — I/I c , and the dimcnsionlcss function L t b is 
d 



+ - 



tan 



/sinh z ^-sinh' v 



sinh 2 ^-sinh 2 



dx 



}.(35) 



The solid curves in Fig. 7 show plots of Ltb vs I/Ic f° r 
D/a = 0.1, 0.3, 1, and 10 for the example of d/a = 0.001. 
All the curves meet at F = 1 or I = I c , where c = 0, 
such that 



L tb {l) 



d 

12a' 



(36) 



Expanding the right-hand side of Eq. (f3"5|) through 
third order in powers of F = I / I c leads to the approxi- 
mation 



Ltb 



12a 



(1.6855aF 2 + f3F 3 ), 



(37) 



where part of the first term, l.OOOOaF 2 , comes from the 
term 1 — c/a, and the other part, 0.6855aF 2 , comes from 
expansion of the integral in Eq. (|35[) and the result 



16 



(tan -1 u) 3 /u 3 du = 0.6855. 



The factor j3 in Eq. |37|) is given by 



8 (na/D) 
it 2 tanh(7ra/Z?) 



(38) 



(39) 



where K(k) and E(k) are complete elliptic integrals of 
the first and second kind of modulus k = tanh(7ra/_D). 
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FIG. 7: Plots of the dimensionless function describing top- 
and-bottom losses Ltb vs F = 7/7 c for d/a — 0.001 and D/a 
= 0.1, 0.3, 1, and 10 (bottom to top). The solid curves were 
obtained from Eq. (|35[) . and the corresponding dashed curves, 
which on this scale are indistinguishable from the solid curves 
except for D/a — 10, were obtained from Eqs. (|15[1 . (|37[) , and 



When D/a » 1, » 2(a/D) 2 , and in the limit D/a -> 0, 
/3 — > 1. The dashed curves in Fig. 7 show plots of L t b 
obtained from Eq. (gT]) for D/a = 0.1, 0.3, 1, and 10 (bot- 
tom to top). The figure shows that this approximation 
is excellent for all F = I/I c except for D/a = 10 close to 
F = 1. In the limit D/a — > 0, when a — ► and /3 — > 1, 
the approximation in Eq. (|3"T|) becomes exact, yielding 



J*6 



12a 



_F 8 



(40) 



which holds for all F. 



C. Comparison of edge and top-and-bottom losses 



FIG. 8: Plots of the dimensionless function L = L e + Ltb 
describing the sum of the edge losses L e [Eq. 1)210 ] and top- 
and-bottom losses Ltb [Eq. (f35)l ] vs F = 7/7 c for 7)/a = 0.1, 
0.3, 1, and 10 (bottom to top) and d/a = 0.001. The solid 
curves display those portions of L for which F > Fx and 
L e > Ltb, and the dashed curves display those portions for 
which F < Fx and L e < L t t- The dotted curves show L 
calculated from the sum of the approximate expressions for 
L e [Eq. (Unj] and L tb [Eq. (37)] . 
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The total energy dissipated per cycle per strip per unit 
length is Q' = Q' e + Q' tb = /j. I^L(F), the sum of the 
edge and top-and-bottom losses. Figure 8 shows plots of 
L = L e + L tb vs F = I/I c for D/a = 0.1, 0.3, 1, and 
10 and d/a = 0.001, the cases considered in Figs. 6 and 
7. It is useful to define Fx as the value of F where the 
edge and top-and-bottom losses are equal, i.e., where the 
curves of L e vs F and L tb vs F cross. Equations (|2"Tj) 
and (HO) yield for d/a = 0.001 the values F x = 0.304, 
0.096, 0.051, and 0.037 for D/a = 0.1, 0.3, 1, and 10. 
The solid portions of the curves in Fig. 8 show where 
F > Fx and L e > L t b, and the dashed portions show 
where F < Fx and L e < L t b- The edge losses are most 
important when D/a is large, and the top-and-bottom 
losses grow in relative importance when D/a is small. 
The dotted curves show L — L e + Ltb calculated from the 
power-law approximations given in Eqs. (|26p and (|37p. 

Figure 9 displays plots of Fx vs D/a for values of d/a = 
0.0001, 0.0003, 0.001, 0.003, and 0.01. These plots show 
that when D/a is small enough and d/a is large enough, 



FIG. 9: Plots of F x , the value of F = I/I c for which the edge 
losses [Eq. (|21[) ] are equal to the top-and-bottom losses [Eq. 
(|35[) ]. vs D/a for several values of d/a. For given values of D/a 
and d/a, when F > Fx, edge losses exceed top-and-bottom 
losses, but when F < Fx, top-and-bottom losses exceed edge 
losses. 



the top-and-bottom losses exceed the edge losses for all 
values of F = 7/7 c . 



IV. DISCUSSION 

In Sec. II of this paper, I derived general expressions 
for the magnetic-field and sheet-current-density distribu- 
tions for an infinite bifilar stack of thin superconducting 
strips of width 2a, thickness d, and separation D, all car- 
rying current of magnitude 7, but with adjacent strips 
carrying current in opposite directions. The calculations 
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assumed that B = fioH and that the critical current 
density J c of each strip was uniform and independent of 
the local magnetic flux density. 

In Sec. Ill, I used critical-state theory to derive ex- 
pressions for Q 1 , the hysteretic ac loss per cycle per unit 
length in each tape, where Q 1 is the sum of edge losses 
Q' e and top-and-bottom losses Q' tb . The top-and-bottom 
losses grow in relative importance as D becomes less than 
a. I expressed Q 1 in terms of integrals of the x and y com- 
ponents of the magnetic fields given in Sec. II. Although 
these integrals can easily be evaluated by numerical in- 
tegration, I also expanded Q' in powers of F = I/I c to 
obtain some useful analytic approximations. 

The behavior of the losses in a bifilar stack are very 
different from those in an infinite stack of strips of sepa- 
ration D when the strips all carry the same current in the 
same direction. In the latter case, Q' , the hysteretic loss 
per unit length per cycle in each tape, increases rapidly 
a,s D/a decreases^ The additive effect of the magnetic 
fields generated by nearby strips greatly increases the 
magnetic field at the edges of a given strip, such that 
the edge losses are greatly magnified. When D/a <C 1 
and F = 1 (I = I c ), the resulting Q'(l) is then larger 
than that of an isolated strip by a factor of 2.71(a/D). 
For the case of a finite stack, a similar enhancement of 
the losses over those in an isolated strip has been noted 
experimentally- and theoretically^^ 

In strong contrast, in a bifilar stack, where adjacent 
strips carry current in opposite directions, the magnetic 
fields generated by nearby strips nearly cancel, and when 
D/a <C 1, the magnetic- field distribution is strongly al- 
tered by the presence of adjacent strips. The perpendic- 
ular component of the field is strongly attenuated, and 
the parallel component becomes nearly uniform across 
the width of the strip. In this case, accounting for both 
the edge and top-and-bottom losses in the stack, Q'(l) 
is smaller than that for an isolated strip by a factor of 
0.552(L>/a) 2 + 0.678(d/a) [see Eqs. (23]), (21, and (3J)]. 

Majoros et ali2! used a finite-element method to calcu- 



late the transport ac losses in finite stacks of supercon- 
ducting tapes of elliptical cross section carrying mutually 
antiparallel currents at the critical value I c . They found 
that when the tapes were closely spaced, the losses were 
less than when the tapes were far apart. The results 
in Sec. Ill confirm this general conclusion not only for 
I = I c but also for all I < I c . 

Although the results in Sec. Ill for Q' were derived for 
an infinitely tall bifilar stack of infinitely long strips, they 
should provide an excellent approximation to the hys- 
teretic loss per cycle per unit length for superconducting 
tapes of finite length in fault-current limiters consisting 
of noninductively wound pancake coils with a large num- 
ber of bifilar windings, so long as the radius r of each 
winding is much larger than the spacing D between ad- 
jacent tapes. 

A key assumption made in this paper is that the mag- 
netic induction in the superconducting film is given by 
B = HqH. This should be a good approximation in 
a high-K superconducting film when the magnitude of 
the self-field H at the edges or surfaces is much larger 
than the lower critical field H c \. However, if this condi- 
tion is not met, then the loss expressions given in Sec. 
Ill will generally overestimate the true losses. For ex- 
ample, for the case of D/a <C 1, when K z w I /2a, if 
H x (x, — e) = K z /2 m I /4a is less than H c i, no vortices 
will be able to penetrate the top or bottom of the strip. In 
this case, the top-and-bottom losses will be zero, rather 
than what was calculated in Sec. Ill B. 
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